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FLOW AND DRAG FORMULAS FOR SIMPLE QUADRICS 

By A. F. Zahm 


PREFACE 

In this text are given the pressure distribution and resistance found by theory and experi- 
ment for simple quadrics fixed in an infinite uniform stream of practically incompressible fluid. 
The experimental values pertain to air and some liquids, especially water; the theoretical refer 
sometimes to perfect, again to viscid fluids. For the cases treated the concordance of theory 
and measurement is so close as to make a resume of results desirable. Incidentally formulas 
for the velocity at all points of the flow field are given, some being new forms for ready use 
derived in a previous paper and given in Tables I, III. A summary is given on page 536. 

The computations and diagrams were made by Mr. F. A. Louden. The present text is a 
slightly revised and extended form of Report No. 312, prepared by the writer for the Bureau of 
Aeronautics in June, 1926, and by it released for publication by the National Advisory Com- 
mittee for Aeronautics. A list of symbols follows the text. 

PRESSURE AND PRESSURE DRAG 

We assume the fluid, of constant density and unaffected by weight or viscosity, to have in 
all the distant field a uniform velocity q 0 parallel to x ; in the near field the resultant velocity q. 
If now the distant pressure is everywhere p 0) and the pressure at any point in the disturbed flovu 
is 7 + the superstream pressure p is given by Bernouilli’s formula, 

p!Pn = i-qVzo\ (l) 

where p n = pqo 2 l2 f called the “stop” or “stagnation” or “nose” pressure. 

At any surface element the superpressure exerts the drag f p dy d z } whose integral over 
any zone 1 of the surface is the zonal pressure drag, 

D^fpdydz. (2) 

Values of p, I) are here derived for various solid forms and compared with those found by 
experiment. 

PRESSURE MEASUREMENTS 

The measured pressures here plotted were obtained from some tests by Mr. R. H. Smith 
and myself in the United States Navy 8 -foot wind tunnel at 40 miles an hour. Very accurate 
models of brass, or faced with brass, had numerous fine perforations, one at the nose, others 
further aft, which could be joined in pairs to a manometer through fine tubing. Thus the 
pressure difference between the nose and each after hole could be observed for any wind speed. 
Then a fine tube with closed tip and static side holes was held along stream at many points 
abreast of the model, to show the difference of pressure there and at the nose. Next the tube 
was thrust right through the model, to find the static pressure before and behind it. The 
method is too well known to require further description. 

THE SPHERE 

Assume as the fixed body a sphere, of radius a, in a uniform stream of inviscid liquid, as 
shown in Table I. Then by that table the flow speeds at points on the axis x, y and on the 
surface are 

(l-a 3 /arTffo, g,= (l+a72^g a , q t = 1.5q, sin 0, (3) 

where 0 is the polar angle. Figure 1 shows plots of these equations. 

1 A zone is a part of the surface bounded by two planes normal to Qo. Usually one plane is assumed tangent to the surface at its upstream end. 
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To graph p/p n in Figure 1, we subtract from the line y~ 1, first g r 2 /g y 2 to show the pressure 
along x; then q 2 /q 2 to portray the surface pressure. A similar procedure gives the superpressure 
in the equatorial plane. 

The little circles show the actual superpressures found with a 2-inch brass sphere in a tunnel 
wind at 40 miles an hour. These agree well with the computed pressures except where or 
near where the flow is naturally turbulent. 

By (3) and (1), on the sphere's surface p/p n ^ 1-2.25 sin 2 #; hence the zonal pressure drag 
f p.2irydy is 

D — ?ra 2 sin 2 0(l — ^ sin 2 0)p n , (4) 

for a nose cap whose polar angle is 0. With increase of 0, as in Figure 2, D/p n increases to a 
maximum .698 a 2 for 0 = 41° — 50' and /> = 0; then decreases to zero for 0 = 70° — 37'; then to its 
minimum — .3927 a 2 for 0~ir/2; then continues aft of the equator symmetrical with its fore 
part. Thus the drag is decidedly upstream on the front half and equally downstream on the 



Fig. 1.— Velocity and pressure along axes and over surface of sphere; graphs indicate theoretical values: 
circles indicate pressures measured at 40 miles per hour in 8-foot wind tunnel, United States Navy 


rear half, having zero resultant. The little crosses, giving D!p n for the measured pressures, 
show that the total pressure drag in air is dowmstream, and fairly large for a body so blunt as 
the sphere. 

Figure 3 depicts the w hole-drag coefficient 2 C d = 2DIv p a 2 q 0 2 , of a sphere, for the manifold 
experimental conditions specified in the diagram, plotted against Reynolds Number R = 2 qtfi/v, 
v being the kinematic viscosity. For 0.2 <R< 200000, the data lie close to the line. 

C D = 2$R~* :i + .48, (5) 

an empirical formula devised by the writer as an approximation. 

For .5<R <2 (5) fairly merges with Oseen’s formula 

C D = 24R~ l +4.5, (6) 

and for R< .2 Stokes' equation (7 D = 24//? is exactly verified. Both these formulas are theoretical. 
Stokes treated only viscous resistance at small scale; Oseen added to Stokes’ drag coefficient, 
24/ R, the term 4.5 due to inertia. 


From the drag D- Co.S, where S is the model’s frontal area, one derives the drug coefficient Ci,— D;p*S. 
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Over an important R range Figure 3 shows C D = .5, giving as the sphere’s whole drag 

D-.5p m S, (7) 

where S = tt a 2 is the frontal area. That is, the sphere’s drag equals half its nose pressure times 
its frontal area. For R<. 2 Stokes’ value, Z) = 6tt ft a q n , has been exactly verified experi- 
mentally, as is well known. 


y 



Fig. 2.— Pressure and pressure-drag on sphere. Graphs indicate 
theoretical values; circles indicate pressure p/p » measured at 
40 miles per hour; crosses indicate pressure-drag D/p,, com- 
puted from measured pressure 


THE ROUND CYLINDER 

Next assume an endless circular cylinder, of radius a, fixed transverse to the stream, as 
indicated in Table I. By that table the flow speed at points on the axes x, y and on the surface is 

qx = (1 - a 2 /x 2 ) q 0 , q y = (1 + a 2 /y 2 ) q 0 , q t = 2 q 0 sin 6, (8) 

where 8 is the polar angle. Plots of (8) are shown in Figure 4. 

Graphs of plp n > made as explained for the sphere, are also given there, together with experi- 
mental values, marked by small circles, for an endless 2-inch cylinder in a tunnel wind at 40 
miles an hour. The agreement is good for points well within the smooth-flow region. 

y 

On the surface pjp n ^ 1 — 4 sin 2 9. The integral 2 f pdy gives, per unit length of 
cylinder, the zonal pressure-drag formula, 

D/pn *» 2 a sin 0- | a sin 3 0. (9) 

42488 — 27 34 
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Fig. 4. — Velocity and pressures along ajtcs ami over surface of endless cylinder; graphs Indicate theoretical values; 
circles indicate pressures measured at 40 miles per hour in 8-foot wind tunnel. United States Navy 



Fio 5 . Pressure and pressure-drag on endless Fig. 6 .-Drag coefficient for an endless cylinder in steady translation through a viscous fluid 

cylinder. Graphs indicate theoretical values; 
circles indicate pressure j>/p« measured at 40 
miles per hour; crosses indicate pressure-drap 
/)/ p„ computed from measured pressure 
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This is 0, 2a/3 (max.), 0, — 2a/3, for 9 = 0°, 30°, 60°, 90° ; and is symmetrical about the equatorial 
plane x = 0. In Figure 5, the little crosses give D/p n for the measured pressures, and show total 
D/p n = 2.33a. 

Figure 6 delineates the drag coefficient C D plotted against R^= 2 aq 0 jv } from Wieselsberger’s 
(Reference 1) wind tunnel tests of nine endless cylinders held transverse to the steady flow. 
The faired line is the graph of 

C d = 9AR~*+ 1.2, (10) 

an empirical equation devised by the present writer. 

For very low values of R , Lamb derives the formula 

Cd = ( 2.002 - \og c R)[l’ (11) 

whose graph in Figure 6 nearly merges with (10) at = . 3. 

For 1 5000 < R< 200000, Figure 6 gives <7 D = 1.2; hence the drag per unit frontal area is 


0 = 1.2 p n , 


( 12 ) 


which is 2.4 times that for the sphere, given by (7). 

THE ELLIPTIC CYLINDER 

An endless elliptic cylinder held transverse to the stream, as shown in Table I, gives for 
points on x, y and on its surface, 


q z =(l-n)q 0 , ff„=(l +m)q 0 , q t = (1 4- b/a)q 0 sin 6, (13) 

where m, n are as in Table I. Amidships q t = (1 + b/a)q 0 = 2q 0 for a = b, as given by (8). Graphs 
of (13) are given in Figure 7. 

To find a', V for plotting (13), assume a ' and with it as radius strike about the focus an 
arc cutting y. The cutting point is distant b' from the origin. Otherwise, b' = V a ,2 — c 2 y where 
c 2 = a 2 —b 2 = const. 

With a/b = 4 one plots p/p n in Figure 7, as explained for the sphere. The circles give the 
experimental plp n for an endless 2-inch by 8-inch strut, at zero pitch and yaw, in a tunnel wind at 
40 miles an hour. The theoretical and measured pressures agree nicely for all points before, 
abreast, and well behind the cylinder. 

Again, sin 2 0 = a 2 ?/ 2 /(& 4 + c 2 ?/ 2 )i if c 2 = a 2 — 5 2 . Hence on the model 

PlPn = 1 - 2l7?0 2 = 1 - ‘ ^ 14) 

y 

This gives the zonal pressure drag, D = 2f pdy, per unit length of cylinder, oi 


D/p„ = 2y-2(a + b) 2 J 


y 2 d y 
V+chf 


= ~4b 


i -f~ b 


y + 2b 2 


(a±b) 2 


tan 


-ycy 

b 2 ' 


( 15 ) 
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whose graph, for cbjb = 4, appears in Figure 8. It rises from 0 at the nose to its maximum where 
p = 0, then falls to its minimum amidships. 




Whatever the value of ajb, the whole pressure on the front half is negative or upstream, 
as for the sphere and round cylinder, and is balanced by the rear drag. For b fixed it decreases 
indefinitely with b/a . 
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The crosses marking actual values of D/p n found in said test show a downstream resultant 
D. In fact, it is one-third the whole measured drag of pressure plus friction, or one-half the 
friction drag. 

For the cylinder held broadside on, b>a and a 2 — b 2 = —r, lienee changing c 2 to — c 2 under 
the integral sign of (15), we find 


D/ P ,=-4b a y y - 


b 2 


(a 


+ &) 2 b 2 r C U 

c 3 b 2 -cy* 


(16) 


where now c 2 = b 2 — a 2 . With 6 fixed, the upstream pressure drag on the front half increases with 
b/a , becoming infinite for a thin Hat plate. It is balanced by a symmetrical drag back of the 
plate. 

Such infinite forces imply infinite pressure change at the edges where, as is well known, the 
velocity can be q = \2p r /p = oo , in a perfect liquid whose reservoir pressure is p T = co . Otherwise 
viewed, the pressure isp, at the plate's center (front and back) and decreases indefinitely toward 
the edges, thus exerting an infinite upstream push on the back and a symmetrical downstream 
push on the front. In natural fluids no such condition can exist. 

THE PROLATE SPHEROID 

A prolate spheroid, fixed as in Table I, gives for points on x, y and the solid surface, respec- 
tively, the flow speeds 

q x ={\-n)q 0t (1 + q t = (1 +fc a ) q 0 sin 0, (16) 



Flo. 9.— Velocity and pressure along axes and over surface of prolate spheroid. Graphs indicate theoretical values; circles 
indicate pressures measured at 40 miles per hour in 8-foot wind tunnel, United States Navy; dots give pressures found 
with an equal model in British test, R, and M. No. 600, British Advisory Committee for Aeronautics 



Fig, 10.— Pressure and pressure-drag on prolate spheroid. Graphs indicate theoretical values; dots indicate 
measured pressure p/p* from Figure 9; crosses indicate pressure-drag Dtp . computed from measured pressure 
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where k a is to be taken from Table II. Graphs of (16) are given in Figure 9, for a model having 
alb = 4, viz., k a = 0.082. 

For this surface p!p n plots as in Figure 10. For a 2 by 8 inch brass model values of p/p n 
are shown by circles for a test at 40 miles an hour in the United States Navy tunnel; by dots for 
a like test in a British tunnel. (Reference 2.) 

By (16), for points on the surface p!p n = 1 -qr!q<?= 1 7 (1 + W 2 sin 2 9. From this, smce 
sin 2 9 = a 2 y 2 j (6 4 -r c 2 y 2 ), the zonal pressure drag f />. 2 tt y d y is found. Thus 


tt« 2 6 4 


b A + c 2 y 2 


D/p n = iry 2 ~- ^2 (1 +fc)*2/*+-y (l+w 2 log. h : 


(17) 


Starting from ?/^0, Djp n increases to its maximum when p — 0, or sin 9 1/(1 +£ a ); then 

diminishes to its minimum for y = b. Figure 10 gives the theoretical and empirical graphs of 
Djp n for a/6 = 4. 

For b fixed the upstream drag on the front half decreases indefinitely with o/a, becoming 

zero for infinite elongation. 

65 OBLATE SPHEROID 


The flow velocity about an oblate spheroid with its polar axis along stream is given by 
formulas in Table I, and plotted in Figure 11, together with computed values of p/p n . No 
determinations of p or D were made for an actual flow. The formula for D/p n is like (17), 
except that c 2 - 6 2 -« 2 , and k a is larger for the oblate spheroid, as seen in Table II. For b fixed 
the upstream drag on the front half increases indefinitely with 6/a. 



The theoretical flow speeds and superpressures for points on the axis of a circular disk fixed 
normal to a uniform stream of inviscid liquid are plotted in Figure 12, without comparative 
data from a test. One notes that the formulas are those for an oblate spheroid with eccen- 
tricity e = l. 

For 1500<tf o a/V <500000, Wieselsberger (Reference 3) finds for the air drag of a thin 
normal disk, of area S f 

D-l.lPnS, m 


or 2.2 times that for a sphere. For aq 0 l» extremely small, theory gives 

D = 5.1 tt ju a qoj (1^) 

ns is well known. Test data for a complete graph, including these extremes, are not yet available. 
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REGIONS OF EQUAL SPEED 

In the flow field g, p are constant where qt 2 ±qn 2 = constant, viz. where 
q 2i q 2 = ^ 1+m ) 2 sin 2 0+ (1 ~n) 2 cos 2 0 = const. 

In particular for the region <l = ([oi this becomes 

_ n 2- n a '* . 
tane= m2 + m = F‘ tan ^ 


( 20 ) 


( 21 ) 



( 22 ) 


Fig. H— Lines of steady flow, lines of constant speed and pressure, for 
infinite frictiODless liquid streaming past a sphere 

For the sphere n as 2m ss fl 3 /r 3 ; hence (21) becomes 

■> 2r 3 — a 3 

UnH ~ 2 2r t -.S^ 

where r = a' = V x 2 + y 2 . The form of this is depicted in Figure 14. 

* tan 0-y/z is the slope of a radial line through the point (x, V) where (21) cuts a confocal curve a’b’, of Table I. knowing a\b' 1 9 to 'locate 
(t, y) draw across the radial line an arc of «'b' by sliding along thex, y axes a straightedge subdivided as in the ellipsograph. The operation rapid 

and easy. 
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A plot of (21) for an elliptic cylinder, fixed as shown in Table I, is given in Figure 16; for a 
prolate spheroid in Figure 17. 

Besides the region (21), having q = q 0) it is useful to know the limit of perceptible disturbance 
say where q 2 fq 0 2 = 1 ± .01. This in (20) gives 

(1 + m) 2 sin 2 <9+ (1 -n) 2 cos 2 0 = 1 ± .01, (24) 

which applies to all the quadrics here studied. Hence 


n 2 — n 0.01 

m 2 + m ± m(2 + m) cos 2 fl’ 


(25) 


A graph of (25) for a round cylinder is shown in Figure 15. Like plots for the other quadrics 



less liquid streaming past a prolate spheroid. Full-line curve q = y„ refers to stream 
parallel to i; dotted curve q = q 4 refers to stream inclined 10° to i 

If in (20) a series of constants be written for the right member, the graphs compose a family 
of lines of equal velocity and pressure, covering the entire flow field. Rotating Figures 14, 17 
about x gives surfaces of q = q 0 . 

COMPARISON OF SPEEDS 

Before all the fixed models the flow speed is q 0 at a great distance and 0 at the nose; abreast 
of them it is q 0 at a distance, and (1 + k a )q 0 amidships. 

The flux of q — q 0 through the equatorial plane obviously must equal q 0 S where S' is the 
body’s frontal area. Hence two bodies having equal equators have the same flux q 0 S, and the 
same average superspeed or average q — q 0 . But the longer one has the lesser midship speed; 
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SYMBOLS USED IN TEXT 


x, y Cartesian coordinates; also axes of same. 

r\ 0 Polar coordinates. 

a __ Angle of attack of uniform stream. 

s. . . . . Length of arc, increasing with 

0 Inclination to x of normal to confocal 

curves in Table I. 

^ Velocity function. 

^ Stream function. 

q Resultant velocity at any point of fluid. 

q 0 Velocity of distant fluid (parallel to x 

axis) . 

q Xt q u Velocity at points on x and y axes (parallel 

to x axis). 

q t Velocity along confocal surface or model 

surface. 

q n Velocity normal to confocal surface. 

p Density of fluid. 

p_ Viscosity. 


v Kinematic viscosity. 

p„__ Nose pressure = p <?« 2 /2. 

p 0 Pressure in distant fluid. 

p Superstream pressure anywhere. 

D Zonal pressure drag = .A p d y d z. 

I ) Whole drag. 

S Frontal area of model. 

C D Drag coefficient= D'pnS. 

R Reynolds number. 

a Radius of sphere, cylinder. 

a, b Semiaxes of ellipse. 

a', b' Semiaxes of confocal ellipse. 

e Eccentricity of ellipse. 

e' Eccentricity of confocal ellipse. 

c Focal distance = ae = a' e' = V° 2 — b 2 

k a Inertia factor (Table II). 

m, n, m,.__ Quantities defined in Tables I, II. 

t Colatitude (see equation 30). 


TABLE I 

Flow functions for simple quadrics 'fixed in a uniform stream of speed g 0 along x positive 


j Symbol defini- 
tions 


Value of functions at any confocal surfaces of semiaxes a', b' 


See diagram A 
(fig. 20) 


See diagram B 
(fig. 20) 


I— ■ 


See diagram C 
(fig. 20) 


Form of quadric 


Sphere 


Circular cylinder 


Elliptic cylinder 


Prolate spheroid 
e = 1 \* a? — b 2 


Velocity function *p j Stream function 4 


Component velocities 

?t» 


— (1 + m) 9oX, where ^ ^ (j — n) fl 0 r/ 2 , where 


*P— — (1 + m) q„x, 


4 = - (l-n) q>y, 


Oblate spheroid 
° 2 


Circular disk 
a = 0, e= 1 


— (H-m)tfpX, 

*=— (1 -n)q 0 y y 

6 a + b 

b a-\-b 

l ~a' a' + b' 

n -b’a' + b' 

-(1 + m) 

+ 

1 

1 

II 

**■ 

log. \t$-v 

1 + e' 2e' 

,og ' l-«' 1-e' 2 

1 + e 2e 

l0g ‘l-€ 1-C 2 

1 + e 2c 2 

log. 

* — (1 + m) q<,x, 


e'b' . , , 

, —sin e 
a 

e a . f 

—sin l e 

i b 

ea . _. 

, —sin 1 e 
b 

n — 

ea . _. 
—sin l e 

! 0 

= — (1 + m)q u x, 

4 = - 2 q ° y 2 ’ 

l (o'~ si " -le ) 

2 / a'b . _. A 

n= -^U rsln € ) 


Differentiation along 
arc s of either figure 
gives: 

d<p dx ... v 
ds =(1+m)q ° 
sin 0, valid for all the ; 
figures; 

..-g 2— <>-•>.. 

cos 6, for the cylin- 
ders; 

<j 0 cos d, for the axial 
surfaces; viz., sphere, 
spheroids, disk. 

! For o', b f = a, b , Table 
II gives m,; whence 
q t — (l + m.) q Q sin 0, 
as the flow velocity on a 
fixed quadric surface. 
q n = ±0 for disk, since 
n= 1 

Remark — both <?, q n 
can be derived from 
either <p or 4* . 

If q u q n =m ax. q i} q n 
on a r b at any other 
point thereof 

qt = q t sin B, <?* = <?.. f 
cos 6 | 


<f, 4, in 


elliptic coordinates, can be found in textbooks; e. g„ « 71,105,108, Lamb’s Hydrodynamics, 4th Ed. 
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hence its outboard speed wanes less rapidly with distance along y. A like relation obtains 
along x from the nose forward. These relations are shown in the velocity graphs of Figures 
18 and 19. A figure similar to 18, including many models, is given in Reference 4. 



Fig. 18.— Supposed graphs of flow speed abreast of endless round and elliptic* 
cylinders of same thickness fixed transverse to an infinite stream of inviscid 
liquid. At great distance flow speed is ?* 



Fig. 19.— Superposed graphs of axial flow sliced before three endless cylinders 1, 2, and 3 (3 osculating 
2), each fixed transverse to an infinite stream of inviscid liquid. At great distance flow speed is?* 


COMPARISON OF PRESSURES 

The foregoing speed relations determine those of the pressures. The nose pressures all 
are p n = pq 0 2 /2; the midship ones are p~p n — (1 + Jc a ) 2 p n . The drag on the front half of the 
model is upstream for all the quadrics here treated; it increases with the flatness, as one proves 
by (15), (17), and is infinite for the normal disk and rectangle. 

APPLICATION OF FORMULAS 

The ready equations here given, aside from their academic interest in predicting natural 
phenomena from pure theory, are found useful in the design of air and water craft. The formula 
for nose pressure long has been used. That for pressure on a prolate spheroid, of form suitable 
for an airship bow, is so trustworthy as to obviate the need for pressure-distribution measure- 
ments on such shapes. The same may be said of the fore part of well-formed torpedoes deeply 
submerged. The computations for stiffening the fore part of airship hulls can be safely based 
on theoretical estimates of the local pressures. The velocity change, well away from the 
model, especially forward of the equatorial plane, can be found more accurately by theory 
than by experiment. The equation (21) of undisturbed speed shows where to place anemometers 
to indicate, with least correction, the relative speed of model and general stream. 

REFERENCES 

1. Wieselsberger, C.: Physicalische Zeitsehrift, vol. 22. 1921. 

2. Joves, R., and Williams, D. H.: The distribution of pressure over the surface of airship model U. 721, 

together with a comparison with the pressure over a spheroid. Brit. Adv. Com. for Aeron. Reports and 
Memoranda No. 600. 1919. 

3. Wieselsberger, C. : Physicalische Zeitsehrift, vol. 23. 1922. 

4. Taylor, D. W.: Speed and Power of Ships, gives a figure similar to 18 but including more models. 1910. 
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TABLE II 

Inertia factors k a * for quadric surfaces in steady translation along axis a in Figure 20 




Prolate spheroid E=ajb 


Elonga- 

Elliptic cylinder, E=a/b 

1 H-« 0 

1 

Oblate spheroid E = bfa 

lion E 

k a = b - 

. eE 2 — E sin l e 


a 

“ , \+e 2e 

IOg * 1— e — 1— e* 

** e- E sin-'e 

1. 00 

1. 000 

0. 500 

0. 500 

1. .50 

. 667 

. 305 

. 803 

2. 00 

. 500 

. 209 

1. 118 

2. 50 

. 400 

. 157 

1. 428 

3. 00 

. 333 

. 121 

1. 742 

4. 00 

. 250 

. 082 

2. 379 

! 5. 00 ! 

. 200 

. 059 

3. 000 

(>. 00 i 

. 167 

. 045 

3. 642 

7. 00 ! 

. 143 

. 036 

4. 279 

8. 00 

. 125 

. 029 

4. 915 

9. (X) 

. Ill 

. 024 

5. 549 

| 10. 00 

. 100 

. 021 

6. 183 

I “ 

. 000 

. 000 

00 


• In this tabje k* = m, of Table I, viz, the value of m when a', b' = o, b. Lamb (FI. and M. No. 623, Rrit. Adv. Com. 
values in the third column above. For mct'on of elliptic cylinder along b axis inertia factor is k»=o/&. 


Aeron.) gives the numerical 
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VELOCITY AND PRESSURE IN OBLIQUE FLOW a 

PRINCIPLE OF VELOCITY COMPOSITION 

A stream q 0 oblique to a model can be resolved in chosen directions into component streams 
each having its individual velocity at any flow point, as in Figure 21. Combining the individuals 
gives their resultant, whence p is found. 


VELOCITY FUNCTION 

Let a uniform infinite stream q 0 of inviscid liquid flowing past a fixed ellipsoid centered 
at the origin have components U, V, W along x, y, z, taken parallel, respectively, to the semi- 
axes, a, b, c; then we find the velocity potential for q 0 as the sum of the potentials <p b , <p t 
for U f V, W. 

In the present notation textbooks prove, for any point (x } y, z) on the confocal ellipsoid 
a' b' c', 

- (1 + m fl ) TJx t (26) 

and give as constant for that surface 


/*oo 

the multiplier of I 
values of <pt, f <p c gives 


m. = abc(l-abc£- af , d ®' % 

being constant for the model, and \ = a' 2 -a 2 . 


(27) 

Adding to (26) analogous 


v? — (1 + m a ) Ux — (1 -fm*) Vy — (1 + m c ) Wz=— (1 4- m)q 0 h ) (28) 

* This brief treatment of oblique flow was added by request after the preceding text was finished. 

* Simple formulas for this integral and the corresponding b, c ones, published by Qreene, R. S. Ed. 1833, are given by Doctor Tuckerman in 
Report No. 210 of the National Advisory Committee for Aeronautics for 1925. Some ready values are listed in Tables III, IV'. 
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where h is the distance of (jc, y, z) from the plane and m a , m b} m c , m are generalized inertia 

coefficients of a ' b' c f for the respective streams U, V , W, q 0 . For the model itself the inertia 
coefficients usually are written k at k b , k c , Jc. The direction cosines of Ji are 


r _ 1 -h m_ a U _ 1+fflj V l + m c ^ 
1 + in qj 1 + th q 0 f X^m^q 0 


(29) 


as appears on dividing (28) by (1 + m)q 0} the resultant of,(l + m a ) U , (1 + m e ) V, (1 + m e ) W. 

EQLT POTENTIALS AND STREAMLINES 


On a f b f c' the plane sections ^ = constant are equipotential ellipses parallel to the major 
section <p = 0, and dwindling fore and aft to mere points, whieh^we call stream poles, where the 
plane (28) is tangent to a' b' c\ If e is the angle between an y^ normal to a' b' c' and the polar 


y 



Fig. 21,— Superposition of streamline velocities for component plane flows parallel 
to axes of elliptic cylinder 


normal, whose direction cosines are L, M, N, we call the line c = const. a line of stream latitude. 
Thus € is the colatitude or obliquity of a surface element of a' b f c' . The line € = 90° is the stream 
equator. This latter marks the contact of a tangent cylinder parallel to the polar normal, viz, 
perpendicular to the plane (28), as in Figure 22. If l , m, n are the direction cosines of any normal 
to a ' b' c' 

cos € — IL + m M + n N. (30) 


Since the streamlines all cut the equipotentials squarely, 3 the polar streamline must run 
continuously normal to the family of confocal ellipsoids a ’ V c f . Hence it forms the intersec- 
tion of a pair of confocal hyperboloids, and at infinity asymptotes a line parallel to q 0 through the 
origin. This straight line may be called the stream axis. Its equation isx:y: z = V :V : IV. 

* On the model, therefore, the streamlines are longitude lines, viz. orthogonals to the latitude lines. 
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COMPONENT VELOCITIES 

At any point of any confocal surface a' b' c f the streamline velocity q , perpendicular to the 
equipotential ellipse there, has components q nt q h respectively, along the surface normal n 
and the tangent s in the plane of q and n . By (28) we have 


2i = 


dtp d h - 


dh ds 


= q t sin c, 


(31) 


where -b<p/bh = (1 + m)<7^</, = max. q t} is the equatorial velocity. By (26) the inward normal 
velocity due to <p a is 

(1 + ».) Vx--m-n.)U, (32) 

o n 

n a being constant on a' b' c’, as may be shown. Similarly, <p*, <p c contribute — m(l — n b )V, — n 
(1 — » C )W; hence the whole normal component is 

g„= -l(l-n a )U-m(l-n b )V-n(l-n c )W=q n cost, (33) 

where q n = [{\-n a YU 2 + (1 - n„) 2 V* + (1 - n c ) J W 2 } 8 = max. q n is the normal velocity at the stream 
poles. Some values of n„ n„ are given in Tables I, III. One also may find (33) as the normal 

derivative of (28). . . 

We now state (28): At any point of a' 1' c' the velocity potential equals q t h, the equatorial 
speed times the distance from the plane of zero potential. Similarly (31) (33) state: At any 
point of a' b' c’ the tangential speed (q t sin t) equals the equatorial speed times the sine of the 
obliquity; the normal speed (q n cos <) equals the polar speed times the cosine of the obliquity. 
This theorem applies to all the confocals, even at the model w r hcre g»~0. 4 

Incidentally the normal flux through a' b’ c' is f q n cos t-dS = </» f dS*, where is the pro- 

jection of S on the plane of 0 = const, and equals the cross section of the tangent cylinder. 
The whole flux through a’ b' c’ is therefore zero, as should be. 

POLAR STREAMLINE 

Some of the foregoing relations are portrayed in Figure 22 for a case of plane flow. Note- 
worthy is the polar streamline or hyperbola. Starting at infinity parallel to g 0> the polar fila- 
ment runs with waning speed normally through the front poles of the successive confocal sur- 
faces; abuts on the model at its front pole, or stop point; spreads round to the rear pole; then 

accelerates downstream symmetric with its upstream part. Its equation q t = 0 = bbjbs can be 

written from (28) 

q t = (1 +m fl ) V sin 6 — (1 + m b ) V cos 0 = 0, or tan ® = ^- r ma jj (34) 


This asymptotes the stream axis y/x= V/U ; for at infinity m a , m b = 0, and tan 9- V/U. Plane- 
flow values of m a , m b are given in Tables I, III. 

All the confocal poles are given by (34); those of the model are at the stops where 


tan 0 = 


1 +k a v ^ a2 V 
i+ic b u b 2 x 


(37) 


Thus on an elliptic cylinder they are where y/z = b*/a ; 3 . V/U; on a thin lamina they are at 
x = ±c cos a, as given in the footnote. Tables II, IV give values of Jc n , Jc b . 

^7 An Analogous theorem obtains also for any other uniform steady stream, say or heat or electricity, that has zero normal component at the bound - 
ary ellipsoid and zero concentration in the flow field. 

* To graph (34) we may use the known relations. , 

0*1/ 0 Lod) 

tan 9 ~b‘*xb' tftna * 


where tan a *=• V! U is the slope of q a or the asymptote lo (34). Thus (34) becomes a'(b • 
tn., rru, reduces to 

1 V — 

c* cos J a c* sin* o. ’ 


(l+Tnk)/(l+m«), which with the tabulated values of 

(36) 


a hyperbola whose semiaies are c cos a. c sin «, c being the focal distance. In this treatment j-o' cos a, jf«b # sin a. a being a fixed eccentric 
angle of the successive confocal elipses. 
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Each angle of attack has its own flow pattern; each its polar streamline given by (34). 
A close-graded family of confocal ellipses and hyperbolas therefore portrays all the poles and 
polar streamlines in the plane ab for all angles of attack. The family can be written 

£ = a'cosar, 2/ = &'sina. (38) 


Thus, giving a', b' a set of fixed values, then a a set, we have the confocal families 


ar 

a ,: 


+ i 

“ Jv/2 


jr 


= \. 


b' 2 c 2 COS 2 a c' sin' a 

the first being ellipses, the second hyperbolas like (36) below. 

Similarly, the locus ^* = 0, or q = qt, is written from (33). With IF=0, 


Its discussion is of minor interest. 


tan 0 = 


\ — UaTJ 
I ~n b V 


DRAG AND MOMENT 


(39) 


(40) 


Formulas for the pressure p all over the simple quadrics here treated are well known, for 
oblique as well as axial flow, and serve to find the drag and moment. For uniform flow the 
resultant drag is zero; its zonal parts can be found as heretofore. The moment about z is the 
surface integral of p(y dy d z — x dz d z), and generally is not zero. 


REGIONS OF EQUAL SPEED ABOUT OBLIQUE MODELS 

Compounding the velocities (31), (33) at any point in the ab plane, as in Figure 22, gives 
for q constant 

q 2 — [(1 + m a ) U sin 0— (l + m*)F cos 0] 2 + [(1 — n 0 ) U cos 0+ (1 — n b )V sin 0] 2 = const. (41) 
In particular for q 2 = TP+V 2 (41) gives 

tan 0 = ~(4±VflC+ A 2 )=“^tan/3. (42) 

where K= V/U , and 

A ~ (1 + Tfi a ) (1 + w &) — (1 — n a) (1 — n P), B = fn a (2 + m a ) — (2 — n^K 2 , Cd = *""^2 ~ ttl + w &). 


y 



Fig. 22.— Polar streamline and component velocities for uniform stream of inviscid liquid about oblique 

elliptic cylinder 


For an elliptic cylinder, as is well known. 

b a + b _ b a + b 

7 ” a_ a' a' + b'’ na ~b'a' + b'’ 


a a + b 
mb ~b / a' + b'’ 


_ a a + b 
Ut> a f a' 4- b' 


which determines A , B , (7, and thence 0 in terms of a' V . Thus, for an endless elliptic cylinder 
of semiaxes a = 4, 6 = 1, yawed 10° to the stream, i. e., F/J7 = tan 10° = .1763, the graph of (42) 
has the form shown full line in Figure 23. This graph takes the dotted form when F = 0, q 0 — U. 


192 



FLOW AND DRAG FORMULAS FOR SIMPLE QUADRICS 

For a prolate spheriod of semiaxes a = 4, 6 = 1, yawed 10°, the graph of (42) is shown in 
Figure 17. 



Fig. 23.— Lines of steady flow, lines of constant speed and pressure, for infinite frictionless liquid 
streaming across endless elliptic cylinder. Dotted curve refers to stream parallel toi; full-lino 
curve 0 = 0# refers to stream Inclined 10° to i 


The two values of tan fi in (42) are 

tan 0, = ^{A + -y/BC+ A 2 ), tan JBC+A 1 ), (43) 

from which are readily derived 

tan tan (0, + &) = ]^rgr ( 44 ) 

(43) give the x-ward inclinations ft, ft, of the asymptotes of the curves g = As can be 
proved, the interasymptote angle ft — ft remains constant as K{— VI U) varies and the asymp- 
totes rotate through J/£(ft-j-ft) about the c axis. 

Thus, with an elliptic cylinder, giving A, B, C their values at eo makes 

tan (ft ft) = 00 , tan (ft + ft) — » (45) 

hence the asymptotes continue rectangular, as in Figure 23, while with varying angle of attack 
they rotate through J^(ft 4- ft). Or more generally one may show that ^ (ft — ft) — 0 .* . ft — ft — 
const. 

A similar treatment applies to the other figures of Table III. For all the cylinders the 
interasymptote angle is 90°; for the spheroids it is 2tan -1 ij2 = 109° — 28' in the ab plane. 
Figure 17 is an example. If the flow past the spheroids is parallel to the be plane the inter- 
asymptote angle for the curves q = q 0 in that plane is obviously unaffected by stream direction. 
It is 90° for infinitely elongated spheroids; 109° — 28' for all others. Excluded from the gen- 
eralizations of this paragraph are the infinitely thin figures, such as disks and rectangles edge- 
wise to the stream, that cause no disturbance of the flow. Passing to three dimensions, we 
note that the asymptotic lines form asymptotic cones having their vertex at the origin. 
42488—27 35 
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SUMMARY 

For an infinite inviscid liquid streaming uniformly, in any direction, past an ellipsoid or 
simple quadric: 

1. The velocity potential at any confocal surface point equals the greatest tangential 
speed along that surface times the distance from the point to the surface’s zero-potential plane. 

2. The tangential flow speed at said surface point equals the greatest tangential speed 
times the sine of the obliquity, or inclination of the local surface element to the equipotential 
plane. 

3. The normal speed at the point equals the greatest normal speed times the cosine of the 
obliquity. 

4. The locus of q = q 0 is a cup-shaped surface asymptoting a double cone with vertex at 
the center. 

5. The vertex angle of this cone is invariant with stream direction; for cylinders it is 
90°, .for spheroids it is 2tan _l V2 = 109° — 28'. 

6. The velocity and pressure distribution are closely the same as for air of the same 
density, except in or near the region of disturbed flow. 

7. The zonal drag is upstream on the fore half; downstream on the rear half; zero on the 
whole. These zones may be bounded by the isobars, e const. 

For the same stream, but with kinematic viscosity v, if the dynamic scale is R = q 0 dj f v, 
d being the model’s diameter: 

8. The drag coefficient of a sphere is 24 /R for 5<. 2; 285“ 85 4- .48 for 0.2<5<200,000 ; 
and 0.5 for 10 4 <C#<O0 5 . 

9. The drag coefficient of an endless round cylinder fixed across stream is 8*75(2.002 — 
log c 5) for 5<.5; approximately 9.4 5 _,8 +1.2 for 0.5<5<200,000; 1.2 for 10 4 <5<200,000. 

10. For 15,000<5<200,000 the drag coefficient of a round cylinder is 2.4 times that for 
a sphere. 
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TABLE III 

Flow functions for simple quadrics in stream V along y positive 
(For all shapes <t>— — (l + ro^) Vy f <?i = (l + wu) V sin e* g„= — (l — ru) V cos «) 


Shape 

nib 

rib 

Sphere 

ej 

« 3 

a' 3 

Circular cylinder 

a 2 
a' 2 

a 2 

a' 2 

Elliptic cylinder 

a a + 6 
6' a' + 6' 

a a-j-b 
a' a' + b' 

Prolate spheroid 

. 1 + e' 2e' 

l0g 'l-«' l-e'> 

. 1 + e' _ , l-2e' 2 

l0g ”l- C ' _2e I-.-.' 

iog.; +e -2. i 1 =^; 
1 — e 1 — e 2 

log. } + e - 2e > + 

1 — e 1 — e 2 

Oblate spheroid 
e=‘v'6*-o» 

c\ 1 — e' 2 — sin~*c' 

“1 + e 2 . , 

c /, , _ strr’e 

V 1 — e 2 

, 1+e' 2 . _. , 

e , -sin ‘e 

y 1 — e 2 f 

1 + e 2 . " 

e , -sin ‘e 

X 1 — e* 


* « is the angle between b' anti any normal to the confocal surface. 

TABLE IV 


Inertia factors for quadric surfaces in steady translation along axis b in Figure 20 


Elonga- 

Ellip. cyl E — alb 
i a 

Prol. spher. E — a/b 

log. ! +e - , 2 \ 

Obi. spher. E — b i 
, e— E sin _l < 

tion E 

kb ~b 

log’+ e -2d- 2 ? 
1—e 1 — e 2 

b= ~eE 2 (e 2 -\-\)-~E \ 

1. 00 

1. 00 

0. 500 

0. 500 

1. 50 

1. 50 

. 621 

. 384 

2. 00 

2. 00 

. 702 

. 310 

2. 50 

2. 50 

. 763 

. 260 

3. 00 

3. 00 

. 803 

. 223 

4. 00 

4. 00 

. 860 

. 174 

5. 00 

5. 00 

. 895 

. 140 

6. 00 

6. 00 

. 918 

. 121 

7. 00 

7. 00 

. 933 

. 105 

8. 00 

8. 00 

. 945 i 

. 092 

9. 00 

9. 00 

. 954 

. 084 

10. 00 

10. 00 

. 960 

. 075 

CO 

CO 

1. 000 

0 


The numerical values in column 3 are given in Lamb’s paper already cited; those in column 4 are given 
substantially by Doctor Bateman, Report No. 163 National Advisory Committee for Aeronautics, 1923. 
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